Nesting, spin-fluctuations, and odd-gap superconductivity in Na-jCoC^ ■ 2/H 2 
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We have calculated the one-electron susceptibility of hydrated Na^CoC^ and find strong nesting 
nearly commensurate with a 2x2 superstructure. The nesting involves about 70% of all electrons 
at the Fermi level and is robust with respect to doping. This nesting creates a tendency to a charge 
density wave compatible with the charge order often seen at x « 0.5, which is usually ascribed to 
electrostatic repulsion of Na ions. In the spin channel, it gives rise to strong spin-fluctuations, which 
should be important for superconductivity. The superconducting state most compatible with this 
nesting structure is an odd-gap triplet s-wave state. 
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Introduction. The recent discovery of superconduc- 
tivity in the layered oxide Nai/3Co02T.4H20 1] is the 
subject of intense experimental and theoretical research, 
despite its relatively low critical temperature, due to a 
combination of properties that suggest the possibility of 
a non-trivial superconducting state and/or a non-trivial 
pairing mechanism. These include unusual magnetic, 
thermodynamic, and transport properties, the apparent 
proximity to both structural and magnetic instabilities, 
and the frustrated triangular Co lattice. @, 0, 0, 

A variety of possible pairing interactions may be rele- 
vant in this system. Structural instabilities have been re- 
ported in the non- hydrated compound [3J, |8( , suggesting 
the possibility of a related soft mode and, correspond- 
ingly, strong electron-phonon coupling. A high polariz- 
ability of water molecules may be responsible for "sand- 
wich" type superconductivity, where the paired elec- 
trons and the pairing bosons are spatially separated^. 
Band structure calculations yield a ferromagnetic ground 
statepl EH and favor an anti-ferromagnetic state over a 
non-magnetic one, presaging long-range spin fluctuations 
that are observed at some dopings [l2j, Finally, and this 
is the central point of the current Letter, nesting proper- 
ties of the calculated Fermi surface suggest the existence 
of strong anti-ferromagnetic spin fluctuations, which are 
bound to play an important role in superconductivity (as 
well as in the normal transport). First and foremost, they 
make a conventional singlet s-wave state rather unlikely. 

Experimental evidence for a pairing symmetry is still 
inconclusive. There are reports of a lack of reduction of 
the Knight shift through the superconducting transition, 
suggesting a triplet state with vector order parameter 
directed out of the plane 0, Q, 0] . Density of states 
probes: [iSR, NMR, and NQR have so far observed a 
non-exponential relaxation rate behavior below T c [lfil 
Il7| . inconsistent with a fully gapped state; some reported 
a Hebel-Slichter maximum near T c [l5|,[lJ|, suggesting a 
coherence peak in the DOS, while others did not find such 
a maximum, possibly because of impurity scattering. 

Here, we calculate the one-electron susceptibility and 
show that it has strong structure in reciprocal space, 
which is not related to crystal symmetry, but is, ac- 
cidentally, nearly commensurate with the lattice. This 



structure is robust with respect to doping and interlayer 
distance, and may even be responsible for the reported 
superstructures (as opposed to an intuitive picture relat- 
ing them solely to Coulomb ordering of Na ions). The 
calculated spin fluctuations appear to be incompatible 
with either singlet or triplet BCS superconductivity, but 
they are fully compatible with so-called odd-gap super- 
conductivity, with the most favorable symmetry being 
triplet s-wave. If this interaction is indeed responsible 
for the superconductivity, then, somewhat similarly to 
MgB2, superconductivity is driven by the Fermi surface 
pockets which have relatively small volume (but large 
density of states). It is the two-dimensionality of the 
electronic structure that makes this possible. 

Electronic structure, nesting and susceptibility. The 
strength of Coulomb correlations in Na 1 / 3 Co02T.4H20 
(subsequently called NCO) is yet unknown, but even in 
metals as strongly correlated as the high-T c cuprates, 
LDA calculations consistently provide accurate Fermi 
surfaces. The one-electron band structure of the par- 
ent compound, Na a Co02, is by now well understood 
[H E3. The Co (i-bands are split by the octahedral 
crystal field of the surrounding oxygens into three t2 g 
bands per layer, well separated from the two e g bands. 
The former further split in the hexagonal symmetry into 
one a\g and two e' g bands. The a\ g and one of the two e' g 
bands cross the Fermi level, forming, respectively, a large 
hexagonal hole pocket around the T point, and six small, 
elliptical hole pockets. The effect of hydration on the 
electronic structure has been shown to be, for all practi- 
cal purposes, related solely to lattice expansion [2fJ, and 
we will therefore use an unhydrated, but expanded com- 
pound for our calculations. The neglected effects of Na 
or H2O disorder would, if anything, further exaggerate 
two-dimensionality. 

We have performed a tight-binding (TB) fit to our 
paramagnetic full-potential LAPW |2lJ band structure 
near E/. The bond length dependence of the TB pa- 
rameters was incorporated as described in Ref.j^], and 
used to analyze Fermi surface dependence on the inter- 
layer distance, c, with the apical O height set to its re- 
laxed value in the hydrated compound. The fits (de- 
tails to be published elsewhere 23]) have an rms error 
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FIG. 1: (color online) The density of states of Nao.3CoC)2 ■ 
1/H2O and of the two bands crossing the Fermi level, in the 
TB model I23T . The ai s band, which carries 2/3 of all holes 
represents only about one third of N(Ef). 

of only 3 mRy in the relevant energy range. This level 
of accuracy and a dense mesh of k-points (over 30,000 
throughout the BZ) were necessary for good resolution 
of the very small FS hole pockets. We use the TB 
Hamiltonian to calculate the one-electron susceptibility 
Xo(q,w) = Xo(q,w) +iXo(q,w), defined as 

Xo(q, uj) = LfOk+q) - /(ek)]/(e k - e k + q -u)-iS), 

k 

where £k is the one-electron energy, / is the Fermi func- 
tion, and the matrix elements are neglected 0. Fig. |3| 
exhibits the nesting structure in %o(q, at uj — > 0, 

important for superconductivity, and also in xd(q, 0). 

With an increased c, as expected, the bands are com- 
pletely two dimensional within the accuracy of our cal- 
culations. The e g > holes get heavier, and comprise ~ 
70% of the DOS at the Fermi energy (Fig. [!}, though 
the total volume of these pockets is half the volume of 
the central a\ g pocket. The latter, which had formed a 
hexagonal prism with moderately flat faces in the parent 
compound, becomes nearly circular and, hence, its con- 
tribution to the susceptibility is practically featureless. 
However, the six elliptical pockets exhibit very good nest- 
ing and because of two dimensionality, their small size 
can only enhance the susceptibility at the nesting vector, 
leaving the deviation from circular cylindrical shape as 
the only factor determining nesting strength. Indeed, if 
they were exactly circular, all three nesting vectors Qi, 
Q2, and Q3, shown in Fig. Olby solid, narrow dashed and 
wide dashed lines, respectively, would nest perfectly. No 
symmetry requirement forces these pockets to be circular, 
nor is there a restriction imposed on their distance from 
the r point. However, due to their small size they are 
nearly perfectly elliptical so that, for example, pockets 
A and D in Fig. nest nearly exactly, and, accidentally, 
the distance between them is just slightly less than half 
of the reciprocal lattice vector G. As a result, a strong 




FIG. 2: The three main nesting types in Nao.3Co02 • y¥LiO. 
The perfect nesting at Qi is shown as a solid line and the two 
imperfect, Q2 and Q3, in narrow and wide dashed lines. 

double-humped peak appears in the calculated Xo (Fig- 
|3J) at all half-integer reciprocal lattice vectors (the dis- 
tance between the humps measures the deviation of the 
nesting vector from G/2). 

Since Qi w G/2, Q 2 and Q 3 are close to G/4, as seen 
in Fig. |3 The corresponding peaks are suppressed by the 
ellipticity of the e g i pockets, which creates misorientation 
between the pockets A and B or A and C. As long as the 
pockets are perfect ellipses, the peaks at Q2 and Q3 have 
the same amplitude. In the real part of x the peak at Qi 
is broadened but still prominent, but the peaks at Q2, 
Q3 are smeared out. 

Ramifications for structure and superconductivity. We 
have now established that there is prominent nesting- 
related structure in both real and imaginary parts of the 
electronic susceptibility at all half-integer reciprocal lat- 
tice coordinates, and weaker, but noticeable structure 
at all quarter-integer coordinates in the imaginary part. 
Because of two-dimensionality [2{|, in the charge chan- 
nel this structure can lead to a Peierls-type charge den- 
sity wave instability, i.e., a structural transition. Indeed, 
various superstructures have been reported, especially at 
x = 1/2, and are often ascribed to charge ordering of Na 
ions. Our results suggest that the C0O2 planes them- 
selves have a tendency toward superstructure formation, 
even without Na ordering. The observed superstructures, 
presumably, are affected by both factors. Finally, there 
are indications 0, that an antiferromagnetic ordering 
may set in parallel to a structural instability, which in 
our picture would correspond to the condensation of a 
spin density wave at a nesting vector. 

At the compositions where the structure in x' does not 
lead to an instability, one may expect soft modes asso- 
ciated with the corresponding wave vectors. However, 
from the point of view of superconductivity, even more 
interesting are the corresponding spin fluctuations, which 
take advantage not only of the structure in the real part 
of xoj but also of the (much sharper) structure of x'o- 
Let us first concentrate on the strongest peak in Xo'; a t 
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humped peaks corresponding to Qi nesting appear along the 
flat edges of the zone boundary, bottom: Xo(q, 0). A temper- 
ature broadening of .001 mRy was used. 



q = Q 1; and consider possible signs of the order param- 
eter on the corresponding pockets (A and D in FigEJ. 
The relevant part of the linearized equation for the order 
parameter can be written as: 

A(k A ,tUJ n ) = T^ ,1/12 A( MD,iOJ n ) (1) 

Here A stands for the order parameter which is scalar 
for singlet, and vector for triplet pairing. The summation 
includes all Matsubara frequencies uj n , the quasiparticlc 
energies are given by £k, and V(q) is the pairing po- 
tential (positive for attraction). For pairing induced by 
phonons, V(q) is always positive, while for spin fluctua- 
tion exchange it is positive for triplet pairing and negative 
for singlet. Since the pockets A and D are related by spa- 
tial inversion, for inversion symmetric order parameters, 
s, d, etc, a solution to Eq. Q exists only if F(Qi) > 0, 
that is, in the triplet channel, and for antisymmetric or- 
der parameters, p etc., only in the singlet channel. Note 
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FIG. 4: The order parameter phase on the relevant sheets of 
the Fermi surface of various odd-gap superconducting states 
in NCO. Solid black lines indicate pairing interactions and 
dashed lines pair-breaking. In the central panel, the phase 
factors are e' s where 5 is indicated in each pocket. 



that such strict selection rules are related to the small 
size of the e' g pockets. Nesting in a large Fermi surface 
that can support line nodes can, in principle, induce a 
singlet d-wsve superconducting state, as discussed, for 
instance, in Ref. [26J. 

The Pauli principle forbids both singlet-p (Sp) and 
triplet-s (Ts) states, as well as Td, unless the order pa- 
rameter is odd with respect to Matsubara time, as dis- 
cussed first by Berezinski [27} for He 3 , and in the con- 
text of solid state by Balatsky and collaborators 28]. 
In other words, spin fluctuations with q = Q x are pair- 
ing for the odd-gap Sp state (oSp) and for the odd-gap 
Ts state (oTs). The relative stability of these must be 
decided by other interactions. Staying within the spin- 
fluctuation scenario, we include now the two other nest- 
ing vectors, Q2 and Q3. Let us first consider the oS 
state. The order parameter is a scalar, and without 
loss of generality we assume it to be constant within 
each of the pockets, with a phase changing from pocket 
to pocket as in Fig The lowest angular momen- 
tum solution allowed for a real order parameter corre- 
sponds to an /-wave, changing among the pockets as 
A A : A B : A c : A D = 1 : -1 : 1 : -1, which im- 
plies several node lines in the induced gap on the pocket 
E. A complex p-wave order parameter is also allowed, so 
that A A : A B : A c : A D = 1 : e in / 3 : e 2 "^ 3 : -1. In this 
case, the induced order parameter on the pocket E varies 
with angle as e 1 ^ , and is nodeless, and thus energetically 
more favorable. Note that as long as the spin-fluctuations 
with vectors Q2 and Q3 couple to electrons with equal 
strength, they do not contribute to pairing at all in ei- 
ther case, canceling each other completely (because of 
the phase factor). 

On the other hand, the oTs state can actually ben- 
efit from the two "minor" nestings. In this case, there 
are no symmetry-implied nodes, and the simplest so- 
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lution has constant gaps over both types of Fermi sur- 
faces. The order parameters are now spinors, which can 
be translated in a standard way into real-space vectors. 
Several solutions are possible with pair spins either in 
or perpendicular to the hexagonal plane, for instance, 
d = const ■ z, where z is the unit vector in c direction, 
or d= const ■ (x+iy). Importantly, for any of the oTs 
states, all spin fluctuations are pairing, including those 
with q = Q 2 and q = Q 3 . Furthermore, an oTs state 
seems compatible with the limited experimental infor- 
mation available on the pairing symmetry. First, such 
a state would not exhibit exponential behavior in DOS- 
probing experiments (specific heat, NMR/NQR, pene- 
tration depth) 0,^3- Second, a coherence peak in the 
NMR relaxation rate T c , if it exists, should be suppressed 
compared to a conventional eSs state. In most experi- 
ments such a peak is not observed, but a rather weak 
coherence peak was seen in Ref. pjlla- Third, specifi- 
cally the d || z triplet s-state, similar to the chiral p-state 
in Sr2RuC>4, implies no change in the in-plane susceptibil- 
ity, as measured by the Knight shift. Finally, 
non-magnetic impurities are not pair-breaking for an oTs 
state, in agreement with observation |29|. 

Conclusion. We calculated the bare susceptibility of 
Na!/ 3 Co02 in the expanded structure, corresponding to 
Na!/ 3 Co02 • J/H2O. We found that the main contribu- 
tion to the density of states comes not from the large 
Fermi surface pocket at the zone center, but from small 
surrounding pockets. These pockets exhibit strong nest- 
ing features, accidentally commensurate with the lattice, 
and suggest proximity to both charge and spin density 
wave instabilities. The structure of the spin fluctuations 
is such that it favors an odd-gap triplet s-wave supercon- 
ductivity, which seems compatible with existing experi- 
ments. The results reported in this Letter are calculated 
for x=0.3, but we have verified that the described nesting 
features remain up to at least x=0.5. 
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